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SUMMARY 


For  nearly  35  years,  the  U.S.  Government  has  been  acquiring  air 
blast,  ground  response,  and  structure  response  data  from  a variety  of  experi- 
ments beginning  with  the  atomic  bomb  tests  in  the  continental  United  States, 
through  the  thermonuclear  events  in  the  South  Pacific  and,  for  the  last 
decade  or  so,  from  high-explosive  events  conducted  in  the  United  States  and 
Canada . 


The  presently  reported  work  describes  the  status  of  an  effort  to 
consolidate  some  of  these  data  in  a repository  known  as  the  DNA  Master  File 
of  Ground-Shock,  Air-Blast,  and  Structure-Response  Data,  but  often  called 
the  DNA  Archive  for  brevity.  A description  of  past  activities  is  presented 
in  this  document,  together  with  the  results  of  a limited-use  program  to 
introduce  the  workings  of  the  system  to  the  technical  community.  Computer 
usage  guidelines  for  the  benefit  of  investigators  wishing  to  estimate  the 
cost  of  processing  data  utilizing  the  capabilities  of  the  data  retrieval  and 
processing  system  that  comprises  part  of  the  DNA  Archive  is  also  presented. 

Much  has  been  written  about  the  DNA  Archive,  and  this  document 
makes  use  of  such  published  information;  it  is  not  intended  to  be  a self- 
sufficient  reference  source.  Moreover,  this  document  covers  a period  when 
most  of  the  time  was  dedicated  to  acquiring  and  filing  data  and  exposing 
the  DNA  Archive  to  potential  users.  Only  a relatively  minor  effort  has  been 
expended  on  using  the  data.  Accordingly,  it  is  likely  that  updated  informa- 
tion will  be  made  available  from  time  to  time.  The  user  should  satisfy 
himself  that  more  recent  information  has  not  been  published  on  the  subject. 
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PREFACE 


The  development  of  the  DNA  Archive,  its  subsequent  exposure  to  the 
technical  community,  and  the  results  and  usage  guidelines  were  all  supported 
by  the  Defense  Nuclear  Agency.  Major  R.  Waters  was  the  Contracting  Officer's 
representative  during  the  early  part  of  the  work.  Me  was  succeeded  by 
Lt.  Col.  D.  Burgess  who,  in  turn,  was  recently  succeeded  by  MaJ . D.  Spangler. 

The  initial  effort  began  in  late  1972  under  Contract  No.  DNA001 “73"C-0058 . 

The  work  continued  under  Contract  No.  DNAOOl -75"C-01 5^  and  was  completed 
under  the  present  contract,  DNAOOl “76-C-Ol 00. 
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attendance  and  participation  of  many  individuals  at  the  DNA  Data  Archive 
Seminar.  Their  participation  is  gratefully  acknowledged.  Special  thanks  are 
extended  to  B.A.  Bolt  of  the  University  of  California,  J.L.  Bratton  of  the 
Air  Force  Weapons  Laboratory,  J.D.  Collins  of  the  J.H.  Wiggins  Company, 

H.F.  Cooper,  Jr.  of  R&D  Associates,  P.F.  Hadala  and  R.E.  Walker  of  the  Waterways 
Experiment  Station,  D.E.  Hudson  of  the  California  Institute  of  Technology, 
and  N.  Lipner  of  the  TRW  Systems  Group  for  their  contributions  to  the  program. 
Finally,  the  interest,  support,  and  involvement  of  M.C.  Atkins,  E.  Sevin, 

D.N.  Burgess,  and  D.R.  Spangler  of  the  Defense  Nuclear  Agency  are  appreciated. 

The  program  manager  at  Agbabian  Associates  was  J.A.  Malthan.  He 
was  assisted  by  M.S.  Agbabian,  R.D.  Ewing,  E.M.  Raney,  and  J.W.  Workman. 

K.T.  Dill,  S.S.  McClellan,  S.D.  Nelson,  and  D.G.  Yates  provided  critical 
support  to  the  program. 
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SECTION  1 


INTRODUCTION 

The  DNA  Data  Archive  project  was  formally  initiated  in  the  fall  of 
1972  as  the  result  of  a long-standing  desire  to  consolidate  ground-shock  and 
air-blast  time  history  data  from  nuclear  and  nuclear-simulation  field  tests 
into  an  archive  for  safekeeping  and  from  which  data  could  be  retrieved  on 
demand.  A major  milestone  in  the  program  was  reached  in  the  fall  of  1975, 
when  the  archiving  of  data  from  the  high-explosive  (HE)  events  listed  in 
Table  1-1  was  completed. 

A description  of  these  test  events  and  the  nature  of  the  data  in 
the  DNA  Archive  is  presented  In  Reference  1. 

The  DNA  Data  Archive  comprises  more  than  a data  retrieval  system. 
It  also  includes  computer  software  providing  data  management  functions,  a 
data  processing  capability  to  perform  a variety  of  time  series  processes  on 
the  data,  and  a data  output  system  for  printed,  taped,  punched,  and  plotted 
formats.  These  capabilities  have  been  reported  extensively  in  References  1, 
2,  and  3- 

In  1975,  the  direction  of  the  program  turned  from  a primarily 
archiving  activity  into  an  information  dissemination  function.  The  purpose 
of  the  latter  was  to  educate  potential  analysts  and  others  interested  in 
these  data  to  the  content  and  use  of  the  Archive.  This  was  achieved  by  the 
presentation  of  a seminar  in  October  1976,  in  which  the  content  of  the 
Archive  was  described  including  the  identification  system  applied  to  the 
data  (Ref.  I),  the  directory  system  that  allows  users  to  locate  particular 
files  (Ref.  ^) , examples  of  data  processing  procedures  that  were  reported  in 
a separate  volume  (Ref.  5),  a summary  of  the  events  whose  data  are  represen- 
ted In  the  Archive  (Ref.  6),  and  a series  of  lectures  (listed  in  Appendix  A) 
describing  contemporary  problems  associated  with  nuclear  weapon  effects. 
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TABLE  1-1.  TIME  HISTORIES  IN  THE  DMA  DATA  ARCHIVE 


Event 

Data  Type 

Ground 

Motion 

Stress  or 
Strain 

Ai r Blast 

Structure 

Response 

DISTANT  PLAIN 

X 

PRAIRIE  FLAT 

X 

X 

DIAL  PACK 

X 

X 

MINE  UNDER 

X 

X 

X 

MINE  ORE 

X 

X 

X 

MINERAL  LODE 

X 

X 

MINERAL  ROCK 

X 

X 

MIXED  COMPANY 

X 

X 

X 

X 

MIDDLE  GUST 

X 

X 

X 

X 

MINE  THROW 

X 

100- lb  HOB 

X 

8-lb  HOB 

X 
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Investigators  interested  in  studying  the  background  and  development 
of  the  DNA  Archive  are  encouraged  to  consult  the  referenced  documents. 

1.1  ORGANIZATION  OF  REPORT 

This  document  is  organized  in  two  basic  parts.  The  first  part 
presents  examples  of  data  processing  that  were  requested  by  users  after  pre- 
sentation of  the  seminar.  The  second  part  includes  guidelines  that  inves- 
tigators may  use  in  estimating  the  computer  and  labor  hours  that  will  be 
required  to  perform  certain  typical  data  processing  operations. 
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SECTION  2 


EXAMPLES  OF  DATA  PROCESSING  REQUESTS 

At  the  completion  of  the  seminar,  attendees  were  invited  to  use 
the  DNA  Archive  at  no  direct  cost  to  themselves.  The  purpose  of  the  offer 
was  to  allow  users  to  familiarize  themselves  with  the  Archive  and  to  provide 
DNA  with  information  on  the  projected  costs  and  future  activities  of  the 
program. 


Four  organizations  responded  to  the  offer:  R&D  Associates;  Physics 
International  Company;  Systems,  Science  and  Software;  and  the  Boeing 
Aerospace  Company.  This  document  presents  the  results  of  all  requests  other 
than  that  of  Systems,  Science  and  Software,  whose  request  was  not  fulfilled 
by  direction  of  DNA. 

In  order  to  provide  maximum  usefulness  to  future  users,  the  proces- 
sing requests  are  presented  as  the  original  request,  the  decisions  related 
to  and  the  results  obtained  from  the  actual  processing,  comments  on  the  output 
by  the  requesters,  and  the  concluding  remarks  referring  to  the  processing 
and  the  interchange  of  information. 

2. 1 REQUEST  FROM  R&D  ASSOCIATES 

When  constructing  shock  spectra  without  resorting  to  direct 
computation,  manual  techniques  are  often  used  to  obtain  an  approximation  of 
shock  spectra.  Thus,  if  an  acceleration  time  history  and  its  integrations 
indicate  peak  acceleration,  velocity,  and  displacement  denoted  by  A,  V,  and  D, 
respectively,  an  approximation  of  the  shock  spectrum  can  be  obtained  by 
constructing  an  envelope  from  the  three  straight  line  segments  on  a tripartite 
graph  where  the  magnitude  of  the  lines  are  2A,  1.5V,  and  D,  respectively. 

The  result  is  typified  in  Figure  2-1. 

The  purpose  of  the  present  task  was  to  determine  how  accurate  this 
approximation  is  for  ground-shock  data  obtained  from  tests  conducted  in  the 
clay-like  materials  of  the  Suffield  test  site  in  Canada. 
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PSEUDOVELOCITY,  LOG  SCALE 


1.5V 

OF  CONSTANT  ACCELERATION 

CONSTANT  DISPLACEMENT 


FREQUENCY,  LOG  SCALE 

FIGURE  2-1.  MANUAL  CONSTRUCTION  OF  SHOCK  SPECTRUM 
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2.1.1  DATA  CONSIDERED 


The  data  in  the  study  included  the  following  files; 


PRAIRIE  FLAT  DIAL  PACK 


2681 

3380 

2649 

3384 

2898 

3362 

2666 

3465 

3453 

3457 

The  reader  should  consult  Reference  4 for  the  identification  of 
these  data. 

2.1.2  PROCESSING  INSTRUCTIONS 

The  instructions  for  processing  these  data  were  quite  specific  and 
occurred  as  negotiated  conversations  via  telephone  or  in  person.  Essentially 
no  written  instructions  were  provided.  The  instructions  are  shown  in 
Table  2-1.  The  data  were  to  be  plotted  as  standard  shock  spectra,  i.e.,  as 
pseudovelocity  plotted  against  frequency. 

2.1.3  PROCESSING  ASSUMPTIONS 

The  reader  can  review  the  input  parameters  associated  with  the 
computation  of  shock  spectra  in  Reference  2 (for  the  SHOXVE  option),  and 
note  that  the  proper  execution  of  the  run  requires  consideration  of  more 
constraints  than  are  provided  by  Table  2-1.  Since  the  maximum  spectra  were 
to  be  computed  only  for  the  forced  era  of  the  time  histories,  there  was  no 
need  to  taper  the  rear  of  the  records  and  the  front  was  not  tapered  in  order 
to  preserve  the  actual  pulse  shape.  Knowing  that  the  spectra  were  to  be 
used  to  corroborate  the  factors  discussed  in  Section  2.1,  the  spectra  plots 
were  exactly  sized  to  fit  standard  tripartite  shock  spectra  graph  paper  so 
that  the  final  plotted  data  on  vellum  could  be  overlaid  on  standard  graph 
paper. 
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File 

2681 

2649 

2898 

2666 

3380 

3384 

3362 

3465 

3453 

3457 


TABLE  2-1.  PROCESSING  INSTRUCTIONS 


Time  Duration  of 

Spectra  Frequency 

Dampi ng 

Measurement , msec 

Spacing 

Ratio 

80 

80 

150 

200 

100 

25  Logari thmi cal ly 
Spaced  Frequenci es 

None 

between 

100 

6 and  8OO  Hz 

100 

200 

250 

300 
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All  data  were  decimated  from  their  various  original  sampling  rates 
to  3000  sps.  This  was  performed  to  reduce  computation  time  in  the  shock 
spectra  processors,  and  it  was  performed  via  the  DECILP  (see  Ref.  2)  proces- 
sing option  in  which  the  default  values  for  the  low-pass  filter  were  selected 
(1500  Hz  cutoff,  6-pole  Butterworth  tangent  type)  with  no  preloading"  provided 
since  the  records  indicated  an  acceleration  that  began  with  essentially  zero 
initial  conditions.  See  Section  2.1.6  for  the  implications  of  data  decimation 

2.1.^  OUTPUT 

The  plotted  output  for  the  processing  request  is  presented  in 
Figures  2-2  through  2-11.  Each  figure  shows  the  original  acceleration  time 
history  and  its  integration.  Displacements  were  not  computed  since  the  very 
low-frequency  characteristics  of  the  shock  spectra  were  not  of  immediate 
interest.  Each  figure  also  shows  the  computed  shock  spectrum.  For  conveni- 
ence, the  approximated  envelope  developed  from  the  factors  shown  in  Figure  2-1 
are  overlaid  as  dashed  lines  on  the  figures. 

2.1.5  REQUESTOR'S  COMMENTS 

The  data  presented  in  Figures  2-2  through  2-11  were  submitted  to  the 
requestor.  Aside  from  technical  questions  concerning  the  algorithm  used  in 
the  shock  spectrum  processor,  no  comments  were  received. 

2.1.6  CONCLUSIONS 

An  examination  of  the  plots  in  Figures  2-2  through  2-11  suggests 
that  the  acceleration  records  contained  some  offset  and  drift.  This  is 
especially  evident  in  Figure  2-4  where  the  maximum  value  of  the  velocity 
time  history  occurs  at  late  time,  apparently  as  the  result  of  a baseline 
error  in  the  acceleration  measurement.  I t i s to  be  noted,  however,  that  the 
shock  spectrum  for  this  record  does  not  differentiate  baseline  errors  from 
true  data,  since  the  constant  velocity  line  on  the  spectrum  estimated  from 


Note:  Preloading  is  implemented  if  the  record  to  be  filtered  begins  with 

nonzero  initial  conditions.  It  "warms  up"  the  filter  and  sets  the 
initial  conditions  to  the  measured  values. 
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RV  1051  PRfllRlE  FLAT  MTS 
DECILP  3.0000E*00  O.OOOOE< 


84-1-flV  HES 

-4.3387E*05  1 .8S59E*0S 

“4.3387E+0S  1 .0S59E»OS 

00  0.0000c«00  O.DOOOE*00  O.OOOOE*00 
0 4924  2681 


CURVE  2 PRGE  3 
39FEflV001R23aR0e4  V LI 
3. 00005^03 

O.OOOOE+00  7.9993E-02 
RV  1051  PRfllRlE  FLAT  MTS 
DECILP  3.0000E«00  O-OOOOE^OO 

PUT  a.ooooc«oo 


84-1-AV  HES 

S.9321E*02  5.8996E*02 

■S.9321E>02  S.8996E*02 

o.ooooE»oo  o.ooooE^oo  o.ooooc«oo 

0 4924  2681 

0 4924  268! 


(a)  Acceleration  and  vel oc i ty--t i me  histories 


FREQUENCY,  HZ 

(b)  Shock  spectrum 


FIGURE  2-2.  PLOTTED  OUTPUT  FOR  FILE  2681 
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VEL.  i IPS)  *\0°  ftCC.  f iPSS)  *10' 


PAGE  1 


300.0 

200.0 

100  .0 

0.(1 

-100.0 

■:200.U 

-300.0 

0-0  10-0  20.0  30.0  40.0  50.0  60.0  70.0  80.0 

TIME  (SEC)  *10'^ 


CURVE  1 PAGE  1 
39FEflVO0lfl23OR140  A L 
3 .U0QaE«a3 

O.UOOOE^OO  7.9999E-02 
•AV  1071  PRAIRIE  FLAT  MT4 
DECILP  3 .OOOUE«0O  O.OOOOEtOO 


140-1-fiV  HES 

'3.0056Et05  3.6037t»0S 

~3.0056E»0S  3.6037E»05 

0.aOOOC>UO  O.UQQOC.OQ  O.UOOOCtOO 
0 4922  2649 


CURVE  2 PAGE  1 
39FEfiV0UlR230R140  V LI 
3.a000E*03 

O.OOOOE*00  7.9999E-02 
-flV  1071  PRAIRIE  FLAT  ttT4 
DEciLP  3.uo(ia£«uo  a.ooooE»oa 

PINT  O.OOOQE^Oa 


140- 1-HV 
-4.8715Et02 
-4.a71SE*02 


3.7098E»02^^ 

3.7098E*02 


o.oaooE*oo 

a 


O.OOOOE.OO  O.UOOOFtUO 
4922  2649 


4322  2649 


(a) 


Acceleration  and  veloci ty--t ime 


hi  stories 


(b)  Shock  spectrum 

FIGURE  2-3.  PLOTTED  OUTPUT  FOR  FILE  26^9 
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V£L-  (IPS)  ■10°  ROC.  (IPSS)  "lO' 


PflQE  4 


TIME  (SEC) 


CURVE  1 PRGE  4 

39FEflV001fl230R330  H L 330-1-flV  HES 

i.0000E»03  -3.8412E*04  2.70I7E*a4 

0.a000E*00  l.SOOOE-O!  -3.84l2Et04  2.7017E*04 

-flV  1 121  PRfilRIE  FLAT  MT2 

OECILP  3.u000£«00  O.UOOOE»Oa  Q.OOOOE^OO  a.U00aC*00  O.UOQOE'UO 

0 492S  2698 


U.O  25.0  50. Q 75.0  100. Cl  )2&.0  150.0 

TIME  (SEC) 


CURVE  2 PAGE  4 
33FEflV001fl230R330  V LI 
3.0000E«03 

O.OQOOE^aO  1 .5000E-OI 
-HV  I !21  PRfilRIE  FLAT  f1T2 
occiLP  3.uoaaE»oa  o.uoooEtOo 

riNT  o.aoooEtOO 


330-1-AV  HES 

“2.&588E*01  9.313BE*CI0 

-2.S588E+01  9.3I38E«0D 

o.uoooE«oo  o.uaooE*aa  a.ooouc«ua 

0 492S  2898 


0 4925  2898 


(a)  Acceleration  and  veloci ty-“time  histories 


o 

o 


\ 

\ 

no" 


FREQUENCY,  HZ 

(b)  Shock  spectrum 


FIGURE  2-4.  PLOTTED  OUTPUT  FOR  FILE  2898 
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VEL.  MPSl  «10°  PCC. 


pflGf;  2 


\l 

1 

\ 

.0  25 

.0  50 

.0  7 

r 

5.0  tO( 

iME  ( S 

J.O  125 

EC)  * 

i.o  1S( 
10‘^ 

1.0  17E 

i.o  200 

CURVE  ! PPGE  2 
j9FEflV0Ulfl23DR400  R L 
i.OOOOE*03 

0 .OOOOEtOO  2 .OOOOE-01 
flV  1 151  PRfllRlE  FLAT  MTl 
OECIIP  3.UUOOE*UO  0-UUOOE*QO 


400-1-flV  HES 

~5.342lE»a4  6.b60iE»O3 

-5.3421E*04  6.b663E«03 

0-UaOOE«UO  O.UOOUE*UO  O.UOOOEtUO 
0 4923  2666 


O-U  25.11  50-0  7&.U  100.0  125.0  150.0  175.0  200.0 


TIME  (SEC)  *10‘^ 


CURVE  2 PAGE  2 
39FEfiV001fl230R400  V LI 
3.0000E*03 

o.ooooEtOO  2.aooaE-o) 

-flV  1151  PRfllRlE  FLAT  MTl 
DEC  UP  3 .OOO0C«0O  0 .00005*00 

PINT  0.00005*00 


400-1-flV  HES 

-9.1435E*OI  5.bl96e*01 

-a.USSE^Ol  5.bI98E*0l 

0.00005*00  0-00005*00  0-00005*00 

0 4923  2666 


0 4923  2666 


(a)  Acceleration  and  velocity — time  histories 


FREQUENCY,  HZ 


(b)  Shock  spectrum 


FIGURE  2-5. 


PLOTTED  OUTPUT  FOR  FILE  2666 
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VEL-  ‘(IP3)  -10°  flCC.  flPSS)  »1Q' 


PAGE  6 


200  .0 


100  .0 


0.0 


~J00  .0 


-200.0 


300.0 

0.0  20.0  40.0  GO.O  60.0  100.0 

TIME  (SEC)  *10“^ 


r 

: I 

CURVE  1 PHGE  6 

44FEflVl .5H215R083  HL  63-1  HES 

2.3998E*03  -2.9801E*05  2.2S92E*0S 

0.00O0E»00  9.9675E-D2  -2.98D1E*0S  2.2592E«0S 

3-1  .5-RV  DIAL  PACK 

OECILP  Q.00a0E*O0  O.OOOOE^OO  0.0000c*00  0.0000E*D0  O.OOOOE*00 

0 4927  3380 


0-0  20.0  40-0  00.0  80.0  100.0 


TIME  (SEC)  *10'^ 


CURVE  2 PAGE  6 

44rERVI .5R2I5R0e3  VLI  83-1. 

2.9998E+03  -6.4479E*02 

O.OOOOEtOO  9.9675E-02  -B.4479E*02 

3-1  .5-nV  DIAL  PACK 

OECILP  B.OOOOE*Oa  O.OOOOE»OD  O.OOOOE*00 


~AV  HES 

5 .3646E*02 
S.364€E»02 

O.OOOOE^OO  O-OOOOE* 
4927  33 


(a) 


Acceleration  and 


velocity-time  histories 


(b)  Shock  spectrum 


FIGURE  2-6.  PLOTTED  OUTPUT  FOR  FILE  3380 
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VEL . ( IPS3  wlO°  flCC . ( IPS3)  «10' 


page:  7 


400  .0 

200.0 

0.0 

■200.0 

400.0 

00  20.0  40-0  60.0  60. 0 100. 0 

TIME  (SEC)  MlO"^ 


CURVE  I PAGE  7 
44FEnVl  .5R215R120  A L 
2.9998E*03 

O.OOOOE+00  9.9675E-02 
20-1  .&-flV  DIAL  PACK 
DCCdP  6 .00006*00  O.OOOOE*00 


120-1  .S-HV  MES 

-4.7393E*05  S.3548E*05 

-4.7393E*0S  S.3548E»05 

O.OOOOE*00  0.00006‘00  O.OOOOE*00 
0 4926  3384 


CURVE  2 PAGE  7 
44FEHV1  .5R21SR120  V LI 
2.99986*03 

O.OOO0E*0O  9.9675E-02 
20-1 .5-RV  DIAL  PACK 
DECILP  8.0000E*00  0.00006*00 

PIHT  O.0000E*00 


120-1. S-AV  HES 

6.6261E*02  3.J283E*02 

■6.62B1E*02  3.1283E*02 

0.0000E*00  00000£*00  O.ODOOE*OQ 

0 4928  3384 

0 4928  3384 


(a) 


Acceleration  and  velocity-time 


hi  stories 


FREQUENCY,  HZ 


(b)  Shock  spectrum 


FIGURE  2-7.  PLOTTED  OUTPUT  FOR  FILE  3384 

21 


VEL.  r IPS)  m\0^  fiCC.  ( IPS3)  n\0 


PAGE  S 


O.u  20>0  40.0  60-0  QO-O  100.0 


TIME  (SEC)  «10"^ 


CURVE  1 PAGE  5 

44FEflVl .5R215R18a  HL  lBO-1  KES 

2.9998E»03  -3.467IEvOS  1 .8209E*O5 

0.0000E*00  9.9B7SE-02  -3.4671E*05  1 .0289E*OS 

00-1  .S-flV  DlRl  PACK 

OECILP  a.0000E*00  O.00aOE*00  O.OOOOE^OO  D.OOODE^OO  O-OODOE^OO 

0 4926  3362 


TIME  (SEC)  *10"^ 


CURVE  2 PAGE  5 
44FERV1  .SA215RiaO  V LI 
2.9998E+a3  -4 

O.OOOOE*00  9.9B75E-02  -4 

80~J  .5-RV  DIAL  PACK 
OECILP  8.000CE«00  O.OOOOE*aO  0 

PINT  o.ooaoE^ao 


18D-1.5-AV  HES 

.9099E*02  l.G865E»02 

.9099E*02  1.6865E*02 

.0000E*00  O.OOOOE»Oa  0.00a0E«00 
0 4926  3362 


0 4926  3362 


(a) 


Acceleration  and 


veloci ty-- 


time  histories 


(b)  Shock  spectrum 

FIGURE  2-8*  PLOTTED  OUTPUT  FOR  FILE  3362 
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VEL . UPS)  «10°  RCC.  ( IPSS)  *10 


PAGE  10 


' \ 

— ' i>  < 

0-0  25.0  50-0  75.0  JOO.O  125-0  150.0  175.0  200.0 

TIME  (SEC)  *10-^ 


CURVE  2 PAGE  10 
44FERV1  .5fl21SR270  V LI 
2.999gE*03 

0 .OOOOE*00  1 .9968E-01 
70-1  .S-RV  OIRL  PACK 
OECILP  B .00005*00  O.OOOOE*! 

PINT  O.OOOOE*00 


270-1. 5-RV  HES 

-3.2412E*02  S.9613E*01 

-3.2412E*02  S.9613E*01 

0.00005*00  0.00005*00  0-00005*00 
0 4931  346S 

0 4931  3465 


CURVE  1 PAGE  10 

44rERVl .bR2l5R27a  AL  270-1.5-RV  HES 

2.9998E*03  -1.41B2E*Q5  3.BS32E*04 

O.0000E*O0  1 .9968E-0I  -!.4182E*05  3.6532E*04 

70-1  .5-flV  DIAL  PACK 

D5CILP  0-00005*00  0.00005*00  0.00005*00  0-00005*00  0-00005*00 

0 4931  3465 


(a) 


Accel erat i on 


and  veloci ty- 


tlme  histories 


FRE(iUENCY,  HZ 

(b)  Shock  spectrum 


FIGURE  2-9 


PLOTTED  OUTPUT  FOR  FILE  3^65 
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VEL-  (IPS)  *10®  flCC.  (IPS3)  «t0' 


PAGE  9 


SO.O  75.0  100.0  125. 0 150.0  175.0  200.0  225-0  250. 0 


TIME  (SEC)  «10‘^ 


CURVE  1 PAGE  9 
i4FEflV]  .SR215R425  A L 
2r9998E*03 

ST0004E-02  2.49G9E-01 

25-1 .S-AV  DIAL  PACK 
DEC1LP  B.0O0OE‘00  0.0000E*00 


42S-1.S-RV  WES 

-1.59S5E^04  7.4956E*03 

-|.59S5E^04  7.4956E*03 

O.OOaOE^OO  O.OOOOE<00  O.OOOOE«OQ 
0 4930  3453 


50.0  75.0  100.0  125.0  150-0  175-0  200-0  225-0  250.0 


TIME  (SEC)  *10“^ 


CURVE  2 PAGE  9 
44FERV1  .5R215R425  V Li 
2 .9396E*03 

5.0004E-a2  2.4g63E-0! 

25-1  .5-AV  DIAL  PACK 
DCCILP  O-ODOOE^OO  O-OGOOE^OO 

PINT  O-OOOOE^OO 


425-1. 5-AV  HES 

-8.0S84E*01  4.in6E*01 

-8.0584E*01  4.1116E‘01 

O.OOOOE^OO  Q.0Oa0E*00  D.UaOOE'OO 
0 4930  3453 

0 4930  3453 


(a) 


Acceleration 


and  velocity-- 


time  histories 


(b)  Shock  spectrum 


FIGURE  2-10.  PLOTTED  OUTPUT  FOR  FILE  3^53 

2k 


01*  (S'Sd!)  'Oatt  ■ qOI*  (Sdli  ”13 A 


PfiGt  8 


100.0  125.0  ISO.O  175.0  200.0  225.0  250.0  275.0  300.0 


TIME  (SEC)  *IQ"^ 


CURVE  1 PAGE  8 
♦ 4FEflVl  .5fl21SR64S  fl  L 

3 .000IE»03 

1.0033E-01  2.9999E-01 

4 645-1  .b-flV  OlflL  PACK 
OCCILP  4.UUQUE«(J0  0.0000E»p0 


64S-1.5~RV  WES 

-5.2S77E+03  2.0257E»03 

-S.2S77E*03  2.02S7E»03 

a.ooooE*uQ  o.oQOOE^uu  o.uaooc«uo 

0 4929  34S7 


CURVE  2 PAGE  8 
44FEAV1 .5R21SR645  V LI 
3.aooiE+a3 

1 .0033E-01  2 .9999E-0I 

4 645-1  .5-RV  DIAL  PACK 
OCCILP  4.uoooc*uo  Q.uooac^oo 

Pur  a.uoouc*uo 


645-L.5-RV  WES 

-2.5985E*01  1 .07!9E*0J 

-2.5985E*01  l.a719E*01 


0 .00006*00 
0 


QOOOOE'OO  0.00006*00 
4929  3457 


0 4929  3457 


(a) 


Accel erat i on 


and  velocity- 


time  histories 


FREQUENCY,  HZ 


(b)  Shock  spectrum 


FIGURE  2-11 


PLOTTED  OUTPUT  FOR  FILE  3^57 
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the  largest  amplitude  of  the  velocity  agrees  very  well  with  the  computed 
spectrum.  Thus,  the  purpose  of  the  processing  was  achieved.  However,  it  is 
to  be  noted  that  the  computed  spectra  in  Figure  2-k  and  all  other  figures 
are  not  the  true  spectra,  since  baseline  errors  distort  the  spectra  in  the 
constant  displacement  and  constant  velocity  regimes. 


An  additional  point  to  be  noted  has  to  do  with  the  accuracy  of 
computed  shock  spectra  as  a function  of  frequency  and  sampling  rate.  If  f 


s 


and  f^  are  the  sampling  rate  and  frequency  of  Interest  in  the  computation 
of  the  shock  spectra,  the  expected  error  in  computing  the  spectra  will  be 


(Ref.  7): 


(2-1) 


Thus,  when  f^  = 800  Hz  and  f^  = 3000  sps,  e - 10^,  which  is  the  expected 
error  at  the  highest  frequency  of  interest  in  the  analysis. 

2.2  REQUEST  FROM  PHYSICS  INTERNATIONAL  COMPANY 

An  area  of  high  uncertainty  in  nuclear  weapon  effects  technology 
is  the  crater- i nduced  ground  motion.  Physics  International  suggested  that 
data  in  the  Archive  might  be  used  to  develop  scaling  rules  by  using  the  fast 
Fourier  transform  (FFT).  A selection  of  measurements  was  considered  in  an 
attempt  to  investigate  the  problem. 

2.2.1  DATA  CONSIDERED 


The  data  in  the  study  included  the  following  files: 


MIDDLE  GUST  IV 


DIAL  PACK 


4290 


3437 

3431 

3445 


The  reader  should  consult  Reference  4 for  the  identification  of 


these  data. 
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2.2.2  PROCESSING  INSTRUCTIONS 

The  processing  instructions  es  received  on  the  Job  Request  Form  are 
shown  in  Figure  2-12.  Although  an  air-blast  record  was  requested  from  event 
DIAL  PACK  by  telephone,  it  was  agreed  that  it  would  be  replaced  by  a record 
from  MIDDLE  GUST  IV,  whose  measurement  duration  more  nearly  matched  the 
DIAL  PACK  velocity  measurements. 

A study  of  Figure  2-12  indicates  a rather  detailed  set  of  instruc- 
tions for  data  to  be  extracted,  processed,  and  displayed.  The  instructions 
included  requirements  to  filter,  linear  detrend  if  necessary,  integrate,  and 
plot  the  resulting  velocities  and  displacements.  In  addition,  it  was 
requested  that  the  FFT  amplitude  be  obtained  for  the  velocity  records. 

A similar  set  of  instructions  was  provided  for  the  processing  of 
the  air-blast  measurement. 

2.2.3  PROCESSING  SETUP  AND  ASSUMPTIONS 

Prior  to  performing  detrends  of  data,  it  is  customary  to  plot  the 
original  records  and  their  integrals  to  establish  the  character  of  the 
offsets  or  drifts  displayed  by  the  data.  These  plotted  data  are  not  part 
of  the  requested  output,  but  they  are  shown  in  Figures  2-13  through  2-15 
for  instructive  purposes.  To  conserve  computer  costs,  the  velocity  records 
were  filtered  (to  600  Hz)  with  a 6-pole  Butterworth,  low-pass  tangent  filter, 
with  a cutoff  at  600  Hz,  and  decimated  to  a sampling  rate  of  1500  sps.  No 
detrending  or  filter  preload  was  included.  The  interested  reader  will  want 
to  compare  the  parameters  of  the  TFILLP  option  in  Reference  A with  the 
printed  output  in  Figures  2-13  through  2-15. 

A review  of  the  air-blast  measurement  (File  ^290)  in  Reference  8 
suggested  that  no  detrending  would  be  required;  hence,  this  measurement  was 
not  initially  plotted. 
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JOB  REQUEST  FORM 


Page_ 


PART  1 — DATA  TO  BE  RETRIEVED 

(Specify  Absolute  File  Numbers.) 


\/aLocny ; 3^37, 3^3/ . 

^THes>e  Ptl2£  U\/(S>  DffTH\ 

■g3'  S' 

;270'  s' 

S' 


l\lPSL^T 

= 2-70  'J 


D/P^L.  Pp<.H 


:0 


< § i -Nh  ^ ^ 

i4t^' 


PART  2— PROCESSES  TO  BE  PERFORMED 

(Specify  Processing  Options,  If  possible,  Include  Code  Words,  Parameter  Values,  and  time  and  frequency  limits,) 


1 .C  • W> 
C «i  O o 

= ^■2  5 


u mo 
— o i/>  - 
a.  < CM  L 


VSi-octT/ss : 

L UM)PA£.:^ pLT&P ^ ^ Pf^^es3/c^  As 

-XP/PSAP  PerSA^o , JTA  sat$a/PAO  A ^ piusa  f pe-rpeMPfm-. 
3.PPOT  US^ciTy- 

4~,  0B7P>/M  FST  Toi^O  H-T:  f pPPT  AMP/ATUPS 

^.OBTAIN  DfsPJ-AcSPn^HT  f pJ-ST. 

/^\Z\7,uhST  : 

AFT^/^  Fium^/N^  A-BOi/B  6BrAfH  ^ P^T  FFT  /^'IPUTIlOB 


PART  3 — FORM  OF  THE  OUTPUT 

(Indicate  plot,  print,  tape  or  punch  output.  If  plots  are  requested,  specify  coordinate  system.) 


Pi-OT  : UNesP-  LihJeA/Z. 

More  : P r^rsu  of  TSN.  (/o)  Plots  aps  pA(^uesTaD. 


3 


PART  4 — HOW  THE  PROCESSED  DATA  WILL  BE  APPLIED 

(Explain  why  you  have  requested  the  processing  above  and  how  you  Intend  to  use  the  results.) 


AppPPS  FAOUjy  OF  US!  NS  APof/us  ANB  TH£ 
FFAOA&tsny  OF  l/S/Nfs-  FFT  70  OBTAIN  ScAsSJOf  Puus  fop, 
opA-repjNS'  - Inpuc&o  SPpllsUP  MotiotJ. 


FIGURE  2-12  PHYSICS  INTERNATIONAL  PROCESSING  REQUEST 
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FIGURE  2-13.  PRELIMINARY  OUTPUT  FOR  FILE  3^37 
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FIGURE  2-14.  PRELIMINARY  OUTPUT  FOR  FILE  3431 
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FIGURE  2-15.  PRELIMINARY  OUTPUT  FOR  FILE 


Based  on  the  results  of  the  information  in  Figures  2-13  through 
2-15,  it  was  elected  not  to  detrend  File  3^37.  File  3^31  was  detrended  by 
removing  the  mean  value  of  the  preground  shock  signal  from  the  entire  record. 
File  3^^5  was  detrended  with  a 1 i near- 1 east-squares  correction  applied  to 
the  ent i re  record. 

The  rationale  of  the  foregoing  was  based  on  the  fact  that  File  3^37 
was  a record  so  near  to  ground  zero  that  it  was  in  the  region  of  large  perma- 
nent deformation,  and  a detrend  procedure  could  not  be  justified  without  an 
understanding  of  the  physics  of  the  phenomenon.  Since  a detailed  analysis  of 
phenomenology  was  beyond  the  scope  of  the  present  work,  the  record  was 
accepted  as  it  was  recorded. 

The  measurement  on  File  3^31  was  sufficiently  close  to  the  detona- 
tion point  that  permanent  displacement  had  occurred.  Thus,  only  the  detectable 
offset  error  was  removed  from  the  data.  Conversely,  File  3^^5  represented  a 
measurement  sufficiently  removed  from  ground  zero  to  preclude  permanent  dis- 
placement and  the  entire  record  was  detrended,  thus  returning  its  displacement 
to  zero  in  late  time.  All  detrends  were  performed  with  the  DETN  option. 

None  of  the  discussion  above  should  be  construed  as  a final  recom- 
mendation for  detrending  these  data;  rather,  it  should  be  viewed  as  one  of 
perhaps  many  rational  procedures  that  could  be  applied. 

The  air-blast  measurement  (File  4290)  had  its  time  base  increased 
by  the  cube  root  scaling  law  so  that  the  100-ton  MIDDLE  GUST  data  could  be 
compared  directly  to  the  500-ton  DIAL  PACK  data.  Thus,  the  time  base  of 
File  4290  was  increased  (with  the  PTIMMU  option)  by  the  factor 

= 1.708. 

Subsequent  to  the  filtering  and  detrending,  FFTs  of  the  velocity 
and  air-blast  records  were  calculated  using  the  FOUR  option.  The  record 
lengths  used  were  3000  points  or  two  seconds.  A cosine  bell  taper  was 
applied  to  the  last  10^  of  the  data  (see  Ref.  4 for  details  of  FFT  operations). 
The  real  and  imaginary  components  of  the  FFT  in  FOUR  option  were  converted  to 
amplitude  and  phase  in  the  CPOL  option. 
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2.2.4  OUTPUT 


The  plotted  output  for  the  three  velocity  measurements  are  presented 
in  Figures  2-16  through  2-l8,  and  for  the  pressure  measurement  in  Figure  2-19. 
They  show  the  fully  processed  records  and  their  integrations  (for  the  velocity 
measurements),  together  with  the  FFT  amplitudes. 

2.2.5  REQUESTOR'S  COMMENTS 

The  data  in  Figures  2-13  through  2-19  were  presented  to  the 
requestor  together  with  printed  output  too  voluminous  for  inclusion  here. 

The  requestor  indicated  that  he  had  intended  to  have  the  FFT  amplitudes  pre- 
sented in  logarithmic  coordinates  and  would  request  such  scales  in  the  future. 

2.2.6  CONCLUSIONS 

The  work  performed  for  this  Archive  user  required  the  most  extensive 
series  of  operations  in  data  processing.  The  reader  is  encouraged  to  study 
the  results  presented  in  Section  2.2.4,  together  with  the  processing 
discussions  presented  in  Section  2.2.3,  prior  to  performing  similar  operations. 

2.3  FIRST  REQUEST  FROM  THE  BOEING  AEROSPACE  COMPANY 

A standard  request  was  received  from  the  Boeing  Aerospace  Company 
for  the  processing  of  air-blast  and  structure  velocity  measurements  to  support 
the  MX  program.  The  processing  requests  consisted  of  filtering,  detrending, 
Fourier  transformation,  and  the  production  of  shock  spectra. 

2.3.1  DATA  CONSIDERED 

The  data  in  the  study  included  the  following  files: 

MIXED  COMPANY 

4325  2486 

4326  2558 

4725  2549 

4728  2463 

The  reader  should  consult  Reference  4 for  the  identification  of 


these  data. 
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F^GURE  2-16.  FINAL  OUTPUT  FOR  FILE  3437 
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FIGURE  2-17.  FINAL  OUTPUT  FOR  FILE  3^31 
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FIGURE  2-18.  FINAL  OUTPUT  FOR  FILE  3hh5 
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FIGURE  2-19.  FINAL  OUTPUT  FOR  FILE  4290 


2.3.2  PROCESSING  INSTRUCTIONS 


The  processing  Instructions  as  received  in  the  Job  Request  Form 
are  shown  in  Figure  2-20.  Due  to  program  constraints,  not  all  of  the  data 
requested  to  be  processed  could  be  included  and  only  those  listed  in 
Section  2.3.1  were  considered. 

The  instructions  for  processing  the  data  are  very  complete.  BAC 
foresaw  the  need  to  detrend  and  requested  plots,  the  generation  of  Fourier 
spectra,  and  in  some  cases,  the  production  of  shock  spectra. 

2.3.3  PROCESSING  SETUP,  ASSUMPTIONS,  AND  OUTPUT 

As  is  the  usual  practice,  all  raw  records  were  integrated  to 
establish  the  characteristics  of  baseline  errors  in  the  data.  Prior  to 
plotting  the  original  data  and  performing  the  integrations,  they  were  fil- 
tered to  2000  Hz  and  decimated  to  a sampling  rate  of  5000  sps  using  the 
TFILLP  option.  A standard  6-po1e  Butterworth,  low-pass  tangent  filter  was 
used  with  no  preload.  The  results,  shown  in  Figures  2-21  and  2-22  for  a 
selection  of  records,  represent  typical  measurement  characteristics 
encountered  in  the  study. 

Based  on  the  data  presented  in  Figures  2-21  and  2-22  and  in 
Appendix  B,  certain  standard  detrending  procedures  were  implemented.  Thus, 
in  Figures  2-23  and  2-2^,  each  record  was  detrended  via  the  DETN  option  by 
fitting  a sloped  line  to  a section  of  the  back  of  each  record,  extending 
the  line  into  the  early-time  regime,  and  subtracting  the  line  from  the 
record. 


In  addition,  the  preshock  arrival  noise  was  removed  by  replacing 
the  data  with  zeros  (PSTC) . For  example,  in  Figure  2-23,  the  printed  output 
associated  with  DETN  indicates  a line  fitted  to  the  data  from  the  time  0.4  sec 
for  3000  data  points  (i.e.,  to  a time  of  1 sec),  which  was  subsequently  sub- 
tracted from  the  entire  record.  The  PSTC  option  indicates  that  data  were 
replaced  by  zeros  beginning  at  time  zero  and  extending  to  0.066  sec. 


Note:  Due  to  the  volume  of  data  involved,  records  in  their  various  processed 

forms  are  presented  in  Appendix  B. 
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yilV  JOB  REQUEST  FORM 

Page  ^ of  ^ 

PART  1— DATA  TO  BE 

(Specify  Absolute  File 

RETRIEVED 

Numbers. ) 

MIXED  COMPANY 

h3ZS 

'(730 

25'(9 

'(705 

^4326 

'(698 

2'(63 

'(707 

‘*725 

'(699 

2628 

'*775 

'(728 

'(7'(6 

2'(66 

'*777 

2'486 

2595 

1*780 

2558 

26'(0 

1*782 

PART  2— PROCESSES  TO  BE  PERFORMED 

(Specify  Processing  Options.  If  possible,  Include  Code  Words,  Parameter  Values,  and  time  and  frequency  limits.) 


Provide:  (a)  Record  drift  and  offset  correction 

(b)  Fast  Fourier  Transform 

(c)  Shock  spectrum 


«/>  ^ 
UJ  CM 
O 

< <n 


u < («  o ® 

o *0  -o  J, 
>4-  z z c 2 

< 3 3 

U (D  Z V 

< «/» 

W 03  O "* 

— u — w 

O.  < CM  UJ 


PART  3 — FORM  OF  THE  OUTPUT 

(indicate  plot,  print,  tape  or  punch  output.  If  plots  are  requested,  specify  coordinate  system.) 


Plots:  Velocity  versus  time 
Pressure  versus  time 


PART  it--HOW  THE  PROCESSED  DATA  WILL  BE  APPLIED 

(Explain  why  you  have  requested  the  processing  above  and  how  you  Intend  to  use  the  results.) 


i These  data  will  assist  in  determining  the  response  of  the  MX  Closure 

, ^ to  a i r b 1 a s t . 


^ ” 


FIGURE  2-20.  THE  BOEING  AEROSPACE  COMPANY  PROCESSING  REQUESTS 
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FIGURE  2-22.  FILE  2^*86  FILTERED,  TFILLP,  TO  2000  HZ 


PRGE  16 


lO 

OJ 

CO 


o o 


(M  — 

O O O C 

cn  + 

UJ  UJ  CO  , 
tt:  J 
O)  CD  lO  5 



Cvj  rsj  { 


I lU  UJ 


(D  ^ 


UJ  LU  lU 
' O O O 
O O O 

I o o o 
o o o 

D O — 


■ + + CO 

-*  I UJ  LU  OD  , 

X o o cn  ^ 

UJ  tl_  O O t ■ 

> UJ  o o UJ  : 

a:  CD  o o I L 

Z3  CO  • • ar.  ^ 

LJ  O LO  CD 


O UJ 

2:  z </? 

o < — 


Q 

(_>  ac 
o < Qc: 

J-  £5D  < 


h~  a 

— UJ 

u.  o 


< Q 
UJ  z 
_J  < 

OC  f/i 

< Q 
UJ  Z 


:5  Q 


o 

o 

X 

to 

Ul 

ac 

O- 


to  •* 

o 

LH  OC 

• o 

O O 
Ul 
O OC 
I- 

X 
o o 


X o 
o X 

CC  Ul 

u_  Qc: 


O 

UJ 

o 

Ul 

a: 

j— , 

UJ 

o 


LA 

CN 

CA 

-cr 


oOl^  I9d^  3yn993ycj 


fA 

CN 

I 

CN 


LU 

cs: 

ZD 

CD 

u. 


k2 


2 


a:  <NJ  c\j 
o o 
+ + 
UJ  UJ 

o o 


O CD  CD  CO 

O CD  oo  OO 

* 

UJ  (N)  CsJ  CM 

O 

O 

o 

o 

o 


o 

ut 

> 

o 

3; 


Qc: 


UJ 

LU 

oa 


CO 

< 


^ o 

^ M 

(DOOCD  COtDCOC 
O + + O UJ  o u 

iiJ  Uuf  rr  2 / 

to  (D  o ^ 

to  (o  or:  g g 

CM  CM  U_  g g 

to  to  00 

• • o • 

CM  CM  O CM  CM 


o to 

O *000 
.^000 
LO  1 


dJ  3 

tr  LO 


iO 


UJ  lu  LU 

000 
000 
000 
I 000 


o 

2 to  O 1 C30  -- 
O O O CD 
O + + CO 
— O UJ  liJ  CD 

ni  (z>  o cr 

UJ  > o O 1 ^ 3. 

> UJ  o o UJ  z:  ^ 

li:  to  o o 1 u.  UJ 

dJ  CO  • • z:  f-  CD 

O D LO  O 


< Z -J 

o — 
o < a: 

h-  CQ  < 

I—  Q ^ 
“ UJ  o 

^00 


< Q O to) 
UJ  Z O Q_t 
— I < LU 


^3 


FIGURE  2-24.  FILTERED  FILE  2486  DETRENDED 


Figures  2-25  and  2-26  present  integrations  of  the  records  shown  in 
Figures  2-21  through  2-24  to  indicate  the  effect  of  detrending  on  the  data. 

Figures  2-27  and  2-28  show  the  fast  Fourier  amplitudes  of  the 
records  presented  in  Figures  2-23  and  2-24.  The  transform  was  calculated 
(four)  for  a record  length  of  2500  points  (0.5  sec),  with  a cosine  taper 
applied  to  the  end  of  the  records.  As  before,  the  real  and  imaginary  com- 
ponents produced  by  FOUR  were  converted  to  amplitude  and  phase  via  CPOL. 

Only  the  amplitudes  are  presented. 

Finally,  shock  spectra  were  computed  for  motion  records  and  an 
example  is  presented  in  Figure  2-29.  This  figure  shows  a different  format 
of  plotted  output  that  is  used  when  it  is  desirable  to  control  the  size  of 
the  plotted  frame.  This  figure  size  matches  the  format  of  standardized 
tripartite  shock  spectra  grids.  The  shock  spectra  are  calculated  using  the 
SHOXVE  option  in  which  40  logarithmically  equally  spaced  frequency  intervals 
were  computed  for  the  forced  era  of  each  record  and  with  no  damping. 

2.3.4  REQUESTOR'S  COMMENTS 

At  the  writing  of  this  document  no  comments  concerning  the  processed 
data  were  received. 

2.3.5  CONCLUSIONS 

The  work  discussed  for  the  Boeing  Aerospace  Company  closely  paral- 
lels the  work  performed  for  Physics  International  and  R&D  Associates. 

Somewhat  different  plotting  formats  have  been  utilized  to  satisfy  the  pre- 
sumed needs  of  the  user. 

2.4  SECOND  REQUEST  FROM  THE  BOEING  AEROSPACE  COMPANY 

A second  Job  Request  Form  was  received  from  the  Boeing  Aerospace 
Company  for  the  processing  of  air-blast  and  particle  acceleration  and  velocity 
measurements  to  be  performed  in  connection  with  the  Deep  Basing  program.  The 
requested  processing  involved  detrending,  Fourier  transformation,  integration, 
and  the  generation  of  shock  spectra. 
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FIGURE  2-23.  SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  FILE  2486  WITH  NO  DAMPING 
AND  NO  RESIDUAL 


2.4.1  DATA  CONSIDERED 


The  data  in  the  study  consisted  of  the  following  files: 


MINERAL  ROCK 


2923 

2919 

^476 

2924 

2920 

4477 

3061 

3059 

4478 

3062 

3060 

4479 

The  reader  should  consult  Reference  4 for  the  identification  of 
these  data. 

2.4.2  PROCESSING  INSTRUCTIONS 

The  processing  instructions  received  on  the  Job  Request  Form  are 
shown  in  Figure  2-30.  Again,  due  to  program  restraints,  not  all  of  the  data 
requested  to  be  retrieved  and  processed  could  be  included.  The  files  that 
were  processed  are  listed  in  Section  2.4.1. 

As  with  the  first  Boeing  request,  the  instructions  for  processing 
were  quite  detailed  in  that  the  need  for  detrending  was  foreseen  along  with 
provisions  being  made  for  plots  and  integrations,  the  calculation  of  Fourier 
transformations,  and  the  generation  of  shock  spectra. 

2.4.3  PROCESSING  SETUP,  ASSUMPTIONS,  AND  OUTPUT 

Prior  to  performing  the  requested  detrend  of  the  data,  the  original 
records  and  their  integrals,  excluding  air-blast  measurements,  were  plotted. 
These  plotted  data  are  not  part  of  the  requested  output,  but  a representative 
sample  is  shown  in  Figures  2-31  through  2-33  for  instructive  purposes.  Before 
plotting  the  original  data  and  generating  the  integrations,  the  records  were 
lowpass  filtered  (TFILLP)  to  2500  Hz,  the  accelerations  and  velocities  were 
decimated  to  a sampling  rate  of  6000  sps  and  the  air-blasts  were  decimated  to 
5000  sps. 

"Note:  Due  to  the  volume  of  data  involved,  the  complete  set  of  records  in 

their  various  processed  forms  are  presented  in  Appendix  C. 
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Jk 


JOB  REQUEST  FORM 


P«ge_ 


PART  1--DATA  TO  BE  RETRIEVED 

(Specify  Absolute  File  Numbers.) 

MINERAL 

ROCK: 

(absolute  file  no.)  MINE  ORE:  (gage 

no. ) 

2921 

3059 

1-35.5AV20 

1-73UV20 

2922 

3060 

1-35.5AH20 

1-73UH20 

2910 

4476 

1-50.5AV20 

1-98AV20 

2911 

4477 

1-50.5AH2O 

1-98AH20 

2923 

4478 

2924 

4479 

Absolute  file  no. 

3061 

Line  4 

Line  6 

3062 

4480  4482  4490 

A487  41*811 

2919 

4489  4485 

4493  4481 

2920 

PART  2--PR0CESSES  TO  BE  PERFORMED 

(Specify  Processing  Optt 

ons.  If  possible,  Include  Code  Words,  Parameter  Values 

, and  time  and  frequency  limits.) 

For  all  motions  (i.e.,  acceleration,  velocity,  displacement)  provide: 

(1)  Offset  and  drift  correction 

(2)  Fast  Fourier  Transform 

(3)  Integrated  displacements 

(4)  Shock  spectra  of  accelerations 

For  ail — hla«;r  royldp* 

(1)  Offset  and  drift  correction 


V oa  o 

— 13  — 

a.  < iM  ui 


PART  3— FORM  OF  THE  OUTPUT 

(Indicate  plot,  print,  tape  or  punch  output.  If  ploti  are  requested,  specify  coordinate  system.) 


Plots:  Acceleration 
Veloci ty 
bi spl acement 
Pressure 


versus  time 


< 

cni 


PART  4— HOW  THE  PROCESSED  DATA  WILL  BE  APPLIED 

(Explain  why  you  have  requested  the  processing  above  and  how  you  Intend  to  use  the  results.) 


These  data  will  assist  in  the  assessment  of  survivability  of  surface 
egress  for  deep-based  facilities. 


ro  LA 

I 00 

LA 


— lO 

— O 
fD  esr 


FIGURE  2-30.  THE  BOEING  AEROSPACE  COMPANY  PROCESSING  REQUESTS 
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DIS.  fIN)  *10"'"  VEL.  (IPS)  *10“  ^ nPSS)  *10' 


r 

CURVE  1 PAGE  1 
43PERV0IDH100R250  A L 
S.999SE^-03 

a.0000E4.00  5.9988E-Q: 

0-RV  S-3  n[N  R0CK-ERST 
TFILLP  O.OOOOE*OD  4.0000E*00 
O.OOOQE^OO 


250-10-RV-E  WES 

-7.0073E*03  2.23S9E+03 

-7.0073E+03  2.23S9E+03 

2.S00aE«Q3  O.OaOQE-^QO  6.0000E«-00 
4962  2923  2923 


2000.0 

3000.0 

4000.0 

5000.0 

5000.0 


ZOO.O  300. 0 400.0  SOO.O 

TIME  (SEC) 


€00.0 


(a) 


Raw  file  2923  filtered,  TFILLP,  to  2500  Hz  and  decimated 
to  6000  samples/second 


TIME  (SEC)  *10"^ 


CURVE  1 PAGE  2 
43FEAV010R100R2S0  V LI 
5.9995E^03 

O.QOOOE+00  5.9988E-01 
0-flV  5-3  niN  R0CK-EflST 
TFILLP  O.OO0OE*0O  4.DOOOE*0O 
O.OO0OE-»OO 
PINT  0.0000E»0Q 


2S0-10-RV-E  WES 

-1.3446E*01  2.5853E+01 

-1.3446E+01  2.5053E*Q1 

2.Sa0QE«03  O.OOOOE*OQ  6.000QE^OQ 
0 2923  2923 


0 2923  2923 


(b) 


Filtered  file  2923  integrated,  PINT, 


to  obta in  velocity 


CURVE  1 PAGE  3 
43FEflVai0R10QR250  □ LI  I 
5.9995E»03 

O.OOOOE^OO  5.9988E-0! 

0-RV  5-3  MIN  R0CK-ERST 
TFILLP  O.aOO0E»O0  4.0000E«00 
O.OOOOE-OO 
PINT  O.OQOOE«00 

PINT  O.O0OOE*0O 


2SO-10-RV-E  WES 

4.0490E-01  9.0046E*00 

4.0490E-0!  9.0046E*00 

2.bOOOE»03  Q.OQQO£*00  G.OOOOE^OO 
0 2923  2923 

0 -2923  2923 

0 2923  2923 


(c) 

FIGURE  2-31 . 


Filtered  file  2923  double  integrated,  PINT  and  PINT,  to 
obtain  displacement 

SAMPLE  OF  RAW  ACCELERATION  MEASUREMENT  AND  RESPECTIVE  INTEGRALS 
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DIS.  tIN)  ^ ^ (IPS)  J«10' 


CURVE  I PAGE  13 
43FEVVG!Ofi![JOR2SO  V L 
6.0002E*03 

0'.OOOOE*OO  S.g383  E-01 
D-UV  13-3  niN  ReCK-EAST 
TFtLLP  0 .OOOOEtOO  2 .OOOQE-00  2 .SOOQE^Oa 
O.aaOOE^-OO  4966 


2S0-10-UV-E  WES 

-5.3019E*00  3.3920E^00 

-S.3019E*Q0  3.3920E+00 


Q-OQQQE.OO  e.tJOOOE-OO 
3061  3061 


(a)  Raw  file  3061  filtered,  TFILLP,  to  2500  Hz  and  decimated 
to  6000  samples/second 


CURVE  I PAGE  U 
i3FEVV01QR!00R250  D UI 
6.0002E^03 

O.OOOOE-00  5 .998IE-01 
0-UV  13-3  MIN  R0CK~ERST 
TFiLLP  O.OOOQE*O0  2. 00006*00  2.S000E*03 
0.00a0E*00  0 

PINT  0.00006*00 


250-10-UV-E  WES 

9.7916E-02  2.20J5E-02 

9.7916E-02  2.201SE-02 


(b)  Filtered  file  306l  integrated,  PINT,  to  obtain  displacement 


FIGURE  2-32.  SAMPLE  OF  RAW  VELOCITY  MEASUREMENT  AND  RESPECTIVE  INTEGRAL 
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FIGURE  2-33.  SAMPLE  OF  RAW  AIR-BLAST  MEASUREMENT  (FILE  4^*76)  FILTERED,  TFILLP,  TO 
2500  HZ  AND  DECIMATED  TO  5000  SAMPLES/SECOND 


Review  of  the  data  presented  in  Figures  2-31  through  2-33  and  In 
Appendix  C leads  to  the  conclusions  that  the  acceleration  measurements  con- 
tained a trend  in  the  form  of  constant  offset,  whereas  the  velocity  and  air- 
blast  measurements  contained  relatively  no  trend  in  that  portion  of  the  record 
containing  the  signal.  Therefore,  each  acceleration,  as  exemplified  in 
Figure  2-34,  was  detrended  (DETN)  by  fitting  a line  with  no  slope  to  the 
section  of  the  record  containing  the  signal,  and  subtracting  the  line  from 
that  section.  For  example,  In  Figure  2-34,  the  print  associated  with  the 
DETN  option  indicates  a line  with  no  slope  fitted  to  the  entire  segment  from 
0.025  to  0.190  sec,  which  was  subsequently  subtracted  from  the  entire  segment. 

All  the  records  were  then  integrated  over  the  segment  containing 
the  signal,  and  the  results  for  those  data  in  Figures  2-31  through  2-33  are 
presented  in  Figures  2-34  through  2-36.  Comparison  of  the  integrations  in 
Figures  2-31b  and  2-31c  with  those  in  Figures  2-34b  and  2-34c  exemplify  the 
effects  of  the  detrend. 

The  Fourier  transform  amplitudes  that  were  calculated  from  the  data 
in  Figures  2-34  through  2-36  are  presented  in  Figures  2-37  through  2-39.  The 
transform  was  calculated  (FOUR)  with  a cosine  taper  applied  to  the  end  of  the 
records.  As  before,  the  real  and  imaginary  components  produced  by  FOUR  were 
converted  to  amplitude  and  phase  via  CPOL.  Only  the  amplitudes  are  presented. 

The  shock  spectra  were  computed  for  the  motion  measurements  using 
the  SHOXVE  option,  the  results  of  which  are  presented  in  Figures  2-40  through 
2-41.  These  figure  sizes  match  the  formats  of  standardized  tripartite  shock 
spectra  grids.  The  shock  spectra  consist  of  values  at  40  logarithmically 
equally  spaced  frequency  Intervals,  computed  for  the  forced  phase  of  each 
record  and  with  no  damping. 

2.4.4  REQUESTOR'S  COMMENTS 

At  the  writing  of  this  document,  no  comments  concerning  the  proc- 
essed data  were  received. 
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FIGURE  2-35.  FILTERED  FILE  306l  WITH  PRESHOCK  ARRIVAL  NOISE  REMOVED  AND 
INTEGRATED  TO  OBTAIN  DISPLACEMENT 
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INTEGRATED  TO  OBTAIN  IMPULSE 
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REAR  10^  OF  RECORD  COSINE  TAPERED 
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REAR  10^  OF  RECORD  COSINE  TAPERED 
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FIGURE  2-40.  SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  FILE  2923  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  FILE  306l  WITH  NO  DAMPING 
AND  NO  RESIDUAL 


2.4.5  CONCLUSIONS 


The  work 
formed  in  response 


discussed  above  is  essentially  the  same  as  the  work  per 
to  the  first  Boeing  Aerospace  Company  request. 


64 


SECTION  3 


GUIDELINES  FOR  ESTIMATING  LABOR  AND  COMPUTER  COSTS 


The  test  data  acquired  over  decades  has  been  accumulated  in  myriad 
forms  and  stored  by  many  different  organizations.  Providing  the  data  in  usable 
and  coherent  form  to  the  DNA  technical  community  involves  collecting,  archiving, 
retrieving,  and  processing  the  data  in  consistent  formats*  While  the  cost  of 
these  operations  is  not  trivial,  the  benefits  realized  by  having  such  data 
readily  available  are  persuasive  in  continuing  such  an  activity.  Therefore, 
the  provision  of  dependable  cost  estimates  for  data  retrieval  and  processing 
are  prerequisites  for  planning  any  data  processing  activity. 


In  this  section,  guidelines  for  cost  estimating  based  primarily  on 
the  archiving  and  processing  described  in  Section  2 are  examined.  The  foil 
i ng  specific  areas  are  discussed: 

a.  Addition  of  edited  data  to  the  DNA  Archive 

b.  Generation  of  a formatted  data  tape  from  the  DNA  Archive 

c.  Processing  of  DNA  Archive  data  through  a typical  sequence  of 
operations  resulting  in  shock  spectra 


All  computer  time  estimates  are  based  on  the  work  performed  by  the 
AA  technical  staff  on  a UNIVAC  1108,  under  the  EXECUTIVE  II  operating  system. 


3.1  CONSIDERATION  FOR  USE  OF  GUIDELINES 

The  computer  estimation  guidelines  presented  in  Section  3.2  are 
based  on  assumptions  regarding  the  number  of  samples  contained  in  a file  and, 
in  some  cases,  sampling  characteristics.  These  assumptions,  discussed  in 
Section  3.2,  should  be  studied  before  using  the  guidelines.  Refer  to  Appendix  D 
for  development  of  generally  applicable  guidelines. 


'More  general  guidelines  to  estimating  costs  for  other  types  of  archiving  and 
processing  are  presented  in  Appendix  D. 
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3.2  ESTIMATION  GUIDELINES 


The  following  is  a presentation  of  computer  and  labor  estimating 
guidelines  for  several  archiving  and  processing  areas. 

3.2.1  ADDITION  OF  EDITED  DATA 

The  following  guideline  is  considered  very  reliable  when  estimating 
the  computer  time  required  to  add  a selected  number  of  measurements,  each  of 
which  contains  from  35,000  to  55,000  samples,  to  the  DNA  Archive: 

Computer  Time,  sec  = 75  x number  of  data  channels  + 90 

Generally,  the  labor  required  to  perform  this  function  can  be 
estimated  as  follows: 

Labor  Hours  - 0.2  x number  of  data  channels  + 7 

Figure  3-1  graphically  displays  the  estimating  guidelines  for 
addition  of  edited  data  to  the  DNA  Archive.  If  the  number  of  data  samples 
exceeds  the  range  specified  above,  it  Is  recommended  that  the  procedures 
discussed  in  Appendix  D be  implemented. 

3.2.2  GENERATION  OF  FORMATTED  DATA  TAPE 

The  following  guideline  is  considered  very  reliable  when  estimating 
the  computer  time  required  to  generate  a formatted  data  tape  containing 
measurements  of  5000  to  7000  data  samples  from  the  DNA  Archive: 

Computer  Time,  sec  = 100  x number  of  data  channels  + 70 

The  labor  required  can  be  estimated  as  follows: 

Labor  Hours  = 0.35  x number  of  data  channels  + 5 
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LABOR  HOURS  TIME,  SEC  X 10 


0 100  200  300 


NUMBER  OF  DATA  CHANNELS 
(a)  Computer 

NOTE:  NUMBER  OF  SAMPLES  FOR  DATA 
CHANNELS  RANGES  FROM  35,000 
TO  55,000 


(b)  Labor 


FIGURE  3-1.  ESTIMATING  GUIDELINES  FOR  ADDITION  OF  EDITED 
DATA  TO  DNA  ARCHIVES 
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Figure  3“2  graphically  displays  the  estimation  guidelines  for 
generation  of  a formatted  data  tape. 

3.2,3  SHOCK  SPECTRA 

The  sequence  of  operations  resulting  in  shock  spectra  consist  of 
the  fol lowing: 


a.  Filter 

b.  Detrend 

c.  Integrate 

d.  Shock  Spectra 

e.  Plot 


The  guideline  presented  here  is  reliable  when  estimating  the 
computer  time  to  calculate  shock  spectra  based  on  the  following  constraints 

a.  Number  of  samples  per  measurement  ranges  from 
40,000  to  50,000 


b.  The  time  period  of  interest  is  10%  of  the 
measured  period 

c.  The  highest  frequency  of  interest  is  10%  of  the 
samp  1 i ng  rate 

The  guidelines  for  estimating  the  computer  time  for  processing  of 
data  through  a typical  sequence  of  operations  resulting  in  shock  spectra  is 
as  foil ows : 


Computer  Time,  sec  = 310  x number  of  data  channels  + 90 
The  amount  of  labor  required  in  this  area  can  be  estimated  as 

fol 1 ows: 


Labor  Hours  = 6.5  x number  of  data  channels  + 15 
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LABOR  HOURS  TIME,  SEC  X 10 


0 50  100 


NUMBER  OF  DATA  CHANNELS 
(a)  Computer 

NOTE:  NUMBER  OF  SAMPLES  PER 
DATA  CHANNEL  RANGES 
FROM  5,000  TO  7,000 


(b)  Labor 


FIGURE  3-2.  ESTIMATING  GUIDELINES  FOR  GENERATION  OF  FORMATTED  DATA  TAPE 
FROM  THE  DNA  ARCHIVE 
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Figure  3“3  graphically  displays  the  estimating  guidelines  for 
processing  of  DNA  Archive  data  through  a typical  sequence  of  operations  that 
will  result  in  shock  spectra. 
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APPENDIX  A 


List  of  Papers  Presented  at  the 
DNA  Data  Archive  Seminar  in  Los  Angeles 
on  13  and  14  October  1976 
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APPENDIX  B 

Processed  Data  from  the  First  Request 
from  the  Boeing  Aerospace  Company 
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FIGURE  B-1.  RAW  AIR-BLAST  RECORD  FILTERED,  TFILLP,  TO  2000  HZ  AND  DECIMATED  TO 
5000  SAMPLES/SECOND 
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REAR  10%  OF  RECORD  COSINE  TAPERED 
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FIGURE  B-6.  FILTERED  AIR-BLAST  RECORD  DETRENDED,  DETN,  WITH  A LINEAR  LEAST  SQUARES 
FIT  TO  DATA  FROM  TIME  BETWEEN  0.4  AND  1.0  SECONDS  AND  LINE  EXTENDED 
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ARRIVAL  NOISE  HAS  BEEN  REMOVED  FROM  RECORD,  PSTC. 
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FIGURE  B-19.  FILTERED  VELOCITY  RECORD  DETRENDED 
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SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  VELOCITY  RECORD  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  VELOCITY  RECORD  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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FIGURE  B-31.  FILTERED  VELOCITY  RECORD  DETRENDED,  DETN,  WITH  A LINEAR  LEAST  SQUARES 
FIT  TO  DATA  FROM  TIME  BETWEEN  0.6  AND  1.0  SECONDS  AND  LINE  EXTENDED 
BACK  INTO  AND  REMOVED  FROM  TIME  0 TO  0.5  SECONDS.  PRESHOCK  ARRIVAL 
NOISE  HAS  BEEN  REMOVED  FROM  RECORD,  PSTC 


pfiGt:  2 


CD 

in 

•M 


05  05 

U5  LD 
<N  Cst 


UJ  UJ 

CO  <o 
CO  CO 
CO  CO 
^ CO  CO 


o o o o o 

OU5  lo  lo  cn 

> 05  05  05  O 


O OO  00 

cn 

M 

cc: 

rvj 

OD  «M  <NJ  CO 
05  o o 

o I 1 z: 

LiJ  IjJ  O 

o o a: 

CvJ  (Si  !j_ 

r-  r^ 

OJ  CM  O 
• • O 

QO  00  cn 
! I Qc: 

to 


o 

o 

o 

o 

o 


cn  o cn  o cs)  o 
o o o 

U UJ  UJ 

o o o 

O O C3 

o o o 

o o — 

CM  <sj  cn 


o 

o 


o 

o o 


o 


o 


(SI  iO 


UJ  ^ 
CD  :x 
cr  • 
a.  o 


a: 

o X 
1 1 
UJ  >■ 
O [ 
O JO 
CO  - 

o — 


o o Q 

o o o 


•lO  -jo 

o 

1 

o 


tu  UJ  UJ 

o o o 
o o o 
o o o 
o o a 


2;  cn  o : o =>  — 
o o CO  cn 
o + + CO 
^ O UJ  UJ  QD  . 

X o o cr 

UJ  > O O I “1  z 

> UJ  CO  O UJ  ^ (_ 

a;  7D  CO  o I u-  UJ 

X CO  ' ' r;  i-  Q 

U)  CO  LO  o 


C_J 

t— 

CO 


o 

a 


O 

< 


Q_ 

CO 


o 


< 
I — 
CO 

o 


Q- 


Q 


< 

oc: 

CD 


B-hO 


pnGt 


qQ\^  ( 9dl j ‘ 13A 
9nina0w  - j.jj 


B-41 


REAR  10%  OF  RECORD  COSINE  TAPERED 
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SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  VELOCITY  RECORD  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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SHOCK  SPECTRA,  SHOXVE , OF  CLEANED-UP  VELOCITY  RECORD  WITH  MO  DAMPING 
AND  NO  RESIDUAL 
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Processed  Data  from  the  Second  Request 
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FIGURE  C-1.  RAW  ACCELERATION  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED 
TO  6000  SAMPLES/SECOND 
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FIGURE  C-2.  FILTERED  ACCELERATION  RECORD  INTEGRATED,  PINT,  TO  OBTAIN  VELOCITY 
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TO  OBTAIN  DISPLACEMENT 
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FIGURE  C-^.  FILTERED  ACCELERATION  RECORD  DETRENDED,  DETN,  WITH  AN  OFFSET  LEAST 
SQUARES  FIT  TO  DATA  FROM  TIME  BETWEEN  0.025  AND  0.190  SECONDS, 
INTEGRATED,  PINT,  TO  OBTAIN  VELOCITY,  AND  DOUBLE  INTEGRATED,  PINT  AND 
PINT,  TO  OBTAIN  DISPLACEMENT.  PRESHOCK  ARRIVAL  NOISE  HAS  BEEN  REMOVED 
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FIGURE  C-5.  FOURIER  TRANSFORM,  FOUR  AND  CPOL,  OF  CLEANED-UP  ACCELERATION  RECORD 
WITH  REAR  10^  OF  RECORD  COSINE  TAPERED 
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FIGURE  C-6.  SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  ACCELERATION  RECORD  WITH  NO 
DAMPING  AND  NO  RESIDUAL 
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FILTERED  ACCELERATION  RECORD  DETRENDED,  DETN,  WITH  AN  OFFSET  LEAST 
SQUARES  FIT  TO  DATA  FROM  TIME  BETWEEN  0.025  AND  0.190  SECONDS, 
INTEGRATED,  PINT,  TO  OBTAIN  VELOCITY,  AND  DOUBLE  INTEGRATED,  PINT 
AND  PINT,  TO  OBTAIN  DISPLACEMENT.  PRESHOCK  ARRIVAL  NOISE  HAS  BEEN 
REMOVED 
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SHOCK  SPECTRA,  SHOXVE , OF  CLEANED-UP  ACCELERATION  RECORD  WITH  NO 
DAMPING  AND  NO  RESIDUAL 
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FIGURE  C-13.  RAW  VELOCITY  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO 
6000  SAMPLES/SECOND 
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FILTERED  VELOCITY  RECORD  INTEGRATED,  PINT,  TO  OBTAIN  DISPLACEMENT 
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FILTERED  VELOCITY  RECORD  WITH  PRESHOCK  ARRIVAL  NOISE  REMOVED  AND 
INTEGRATED  TO  OBTAIN  DISPLACEMENT 
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SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  VELOCITY  RECORD  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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RAW  VELOCITY  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO 
6000  SAMPLES/SEC 
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FIGURE  C-20.  FILTERED  VELOCITY  RECORD  WITH  PRESHOCK  ARRIVAL  NOISE  REMOVED  AND 
INTEGRATED  TO  OBTAIN  DISPLACEMENT 
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SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  VELOCITY  RECORD  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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FIGURE  C-23.  RAW  ACCELERATION  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO 
6000  SAMPLES/SEC 
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FIGURE  C-25.  FILTERED  ACCELERATION  RECORD  DOUBLE  INTEGRATED,  PINT  AND  PINT,  TO 
OBTAIN  DISPLACEMENT 
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FIGURE  C-26.  FILTERED  ACCELERATION  RECORD  DETRENDED,  DETN,  WITH  AN  OFFSET  LEAST 
SQUARES  FIT  TO  DATA  FROM  TIME  BETWEEN  0.025  AND  0.190  SECONDS, 
INTEGRATED,  PINT,  TO  OBTAIN  VELOCITY,  AND  DOUBLE  INTEGRATED,  PINT 
AND  PINT,  TO  OBTAIN  DISPLACEMENT.  PRESHOCK  ARRIVAL  NOISE  HAS  BEEN 
REMOVED 
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FIGURE  C-28.  SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  ACCELERATION  RECORD  WITH  NO 
DAMPING  AND  NO  RESIDUAL 
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FILTERED  ACCELERATION  RECORD  INTEGRATED,  PINT.  TO  OBTAIN  VELOCITY 
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FILTERED  ACCELERATION  RECORD  DETRENDED,  DETN , WITH  AN  OFFSET  LEAST 
SQUARES  FIT  TO  DATA  FROM  TIME  BETWEEN  0.025  AND  0.190  SECONDS, 
INTEGRATED,  PINT,  TO  OBTAIN  VELOCITY,  AND  DOUBLE  INTEGRATED,  PINT 
AND  PINT,  TO  OBTAIN  DISPLACEMENT.  PRESHOCK  ARRIVAL  NOISE  HAS  BEEN 
REMOVED 
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FIGURE  C-3^.  SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  ACCELERATION  RECORD  WITH  NO 
DAMPING  AND  NO  RESIDUAL 
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RAW  VELOCITY  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO 
6000  SAMPLES/SEC 
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FIGURE  C-37.  FILTERED  VELOCITY  RECORD  WITH  PRESHOCK  ARRIVAL  NOISE  REMOVED  AND 
INTEGRATED  TO  OBTAIN  DISPLACEMENT 
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FIGURE  C-39.  SHOCK  SPECTRA,  SHOXVE,  OF  CLEANED-UP  VELOCITY  RECORD  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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FIGURE  C-40.  RAW  VELOCITY  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO 
6000  SAMPLES/SEC 
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FIGURE  C-k2.  FILTERED  VELOCITY  RECORD  WITH  PRESHOCK  ARRIVAL  NOISE  REMOVED  AND 
INTEGRATED  TO  OBTAIN  DISPLACEMENT 
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REAR  10^  OF  RECORD  COSINE  TAPERED 
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FIGURE  C-^^.  SHOCK  SPECTRA.  SHOXVE,  OF  CLEANED-UP  VELOCITY  RECORD  WITH  NO  DAMPING 
AND  NO  RESIDUAL 
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FIGURE  C-h5.  RAW  AIR-BLAST  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO  5000 
SAMPLES/SECOND 
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INTEGRATED  TO  OBTAIN  IMPULSE 
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REAR  10%  OF  RECORD  COSINE  TAPERED 
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FIGURE  C-48.  RAW  AIR-BLAST  RECORD  FILTERED,  TFI LLP,  TO  2500  HZ  AND  DECIMATED  TO 
5000  SAMPLES/SEC 
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INTEGRATED  TO  OBTAIN  IMPULSE 
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FIGURE  C-50.  FOURIER  TRANSFORM,  FOUR  AND  CPOL,  OF  CLEANED-UP  AIR-BLAST  RECORD  WITH 
REAR  10^  OF  RECORD  COSINE  TAPERED 
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FIGURE  C-51.  RAW  AIR-BLAST  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO 
5000  SAMPLES/SEC 
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INTEGRATED  TO  OBTAIN  IMPULSE 
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FIGURE  C-53.  FOURIER  TRANSFORM,  FOUR  AND  CPOL,  OF  CLEANED-UP  AIR-BLAST  RECORD  WITH 
REAR  10^  OF  RECORD  COSINE  TAPERED 
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FIGURE  C-5^.  RAW  AIR-BLAST  RECORD  FILTERED,  TFILLP,  TO  2500  HZ  AND  DECIMATED  TO 
5000  SAMPLES/SEC 
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INTEGRATED  TO  OBTAIN  IMPULSE 
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FIGURE  C-56.  FOURIER  TRANSFORM,  FOUR  AND  CPOL,  OF  CLEANED-UP  AIR-BLAST  RECORD  WITH 
REAR  10^  OF  RECORD  COSINE  TAPERED 
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General  Labor  and  Computer  Estimation  Guidelines 
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APPENDIX  D 


GENERAL  LABOR  AND  COMPUTER  ESTIMATION  GUIDELINES 

It  has  been  established  that  potential  Archive  users  should  have 
access  to  guidelines  that  will  allow  them  to  estimate  the  computer  hours  and 
labor  hours  required  for  data  processing.  Therefore,  this  appendix  presents 
general  guidelines  resulting  from  an  analysis  based  primarily  on  the  processing 
described  in  Section  2. 

The  guidelines  are  oriented  to  typical  processing  areas 

including: 

a.  Filtering  and  decimation 

b.  Detrending  of  acceleration,  velocity,  and  air-blast  measurements 

c.  Calculation  of  shock  spectra 

d.  Calculation  of  Fourier  transforms 

e.  Calculation  of  both  Fourier  transforms  and  shock  spectra 

All  computer  time  estimates  are  based  on  the  work  performed  with  a 
UN  I VAC  1108,  Executive  I I . 

D.l  REQUIRED  INFORMATION 

In  order  to  perform  an  accurate  estimate,  certain  information  must 
be  compiled.  The  sources  of  this  information  are  to  be  derived  from: 

a.  The  requirements  of  the  user 

b.  Reference  k 

The  following  is  a description  of  the  information  needed  from  each 

source. 
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D.1,1  INVESTIGATOR  PROVIDED  INFORMATION 


The  investigator  must  decide  which  of  the  available  measurements 
are  to  be  processed.  The  following  information  must  then  be  compiled  for 
each  measurement: 

a.  Measurement  type;  i.e.,  air  blast,  velocity,  acceleration,  etc. 

b.  Frequency  limits,  usually  consisting  of  some  lower  limit 

approaching  zero  and  an  upper  limit  representing  the  highest 
frequency  of  interest  (Hz)  , f^ 

c.  Time  limits,  consisting  of  the  length  of  the  time  segment 
of  interest  (sec),  T 

All  the  measurements  should  then  be  combined  into  groups  having 
the  same  measurement  type,  highest  frequency  of  interest,  and  length  of  time 
segment  of  Interest. 

The  information  that  should  be  compiled  for  each  group  consists  of: 

a.  The  number  of  measurements. 

b.  Measurement  type. 

c.  Highest  frequency  of  interest  (Hz),  f . 

m 

d.  Length  of  the  time  segment  of  interest  (sec).  T . 

’ m 

e.  A list  of  the  absolute  file  numbers  that  correspond  to  the 
measurements.  This  can  be  obtained  by  referring  to 
Reference  A or  Reference  6, 
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D.1.2  REFERENCE  PROVIDED  INFORMATION 


Making  use  of  Reference  4,  combined  with  the  list  of  absolute 
file  numbers,  the  user  must  then  compile  information  for  each  measurement 
of  a group  consisting  of  the  following: 

a.  DELTA  X,  the  sampling  increment  (sec),  Ax 

b.  FOOTAGE,  the  location  of  the  measurement  on  tape^  F 

c.  The  previous  entry's  FOOTAGE,  corresponding  to  the  location 
of  the  previous  file  on  the  tape,  F' 

The  combination  of  this  information,  along  with  that  discussed  in 
Section  D.Kl,  is  needed  to  define  the  values  of  the  parameters  used  in  the 
guideline  algorithms. 

An  example  of  an  information  tabulation  for  the  data  discussed  in 
Section  2.1  is  presented  in  Table  D-1. 

D.2  PARAMETER  VALUE  CALCULATIONS 

In  order  to  make  use  of  the  guideline  algorithms  presented  in 
Section  D.3,  the  following  parameter  values  must  be  calculated  for  each 
measurement  of  a group: 

a.  Decimation  ratio,  d 

b.  Number  of  samples,  n^ 

c.  Segment  ratio,  s 

All  groups  are  then  further  subdivided  into  classes  having  the 
same  decimation  ratio,  number  of  samples,  and  segment  ratio  values.  The 
following  is  a discussion  of  the  method  used  to  determine  the  value  of  these 
parameters.  Refer  to  Table  D-2  for  an  example. 
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TABLE  D-2.  SUMMARY  OF  PARAMETER  VALUE  CALCULATIONS  FOR  RDA  REQUEST 


Group 

N umbe  r 
of 

Classes 

Class 

Number 

of 

Measurements , 
m 

c 

Decimation 

Ratio, 

d 

Number 

of 

Samples, 

n 

s 

Segment 

Ratio, 

s 

Abso 1 ute 
File 
Number 

1 

1 

1A 

2 

0.5 

29628 

0.02A 

2681 

2649 

2 

1 

2A 

1 

0.5 

29628 

0.046 

2898 

3 

2 

3A 

1 

0.5 

31668 

0.057 

2666 

3B 

1 

0.2 

96067 

0.050 

3465 

k 

1 

AA 

3 

0.2 

96067 

0.025 

3380 

3384 

3362 

5 

1 

5A 

1 

0.2 

105101 

0.057 

3453 

6 

1 

6A 

1 

0.5 

A8132 

0.075 

3457 
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D.2J  DECIMATION  RATIO 


For  each  measurement,  the  decimation  ratio  is  defined  as  roughly  six 
times  the  sampling  increment  times  the  highest  frequency  of  interest.  This 
can  be  represented  as: 

d = 6(Ax)  f 

m 

where 

d = Decimation  ratio 

Ax  = Sampling  increment  (sec) 

f^  = Highest  frequency  of  interest  (Hz) 

Refer  to  Figure  D-1  for  specific  values  of  decimation  rates. 


D.2.2  NUMBER  OF  SAMPLES 

For  each  measurement,  the  number  of  samples  is  equivalent  to  one 
thousand  four  hundred  fifty-seven  (1,457)  times  the  difference  between  the 
location  of  the  measurement  on  tape  and  the  location  of  the  previous  master- 
file  entry  on  tape  minus  one  thousand  eight  hundred  forty-three  (1,843). 

This  can  be  represented  as: 


= 1457(F  - F‘)  - 1843 

whe  re 

= Number  of  samples 

F = Location  of  the  measurement  on  tape 

F'  = Location  of  the  previous  masterfile  entry  on  tape 


Refer  to  Figure  D-2  for  specific  values  of  number  of 


samples. 


D.2,3  SEGMENT  RATIO 

For  each  measurement,  the  segment  ratio  is  equivalent  to  the  ratio 
of  the  length  of  the  time  segment  of  interest  to  one  less  than  the  number  of 
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DECIMATION  RATIO 


SAMPLING  INCREMENT,  Ax,  SEC 


FIGURE  D-1.  DECIMATION  RATIO  (d)  VERSUS  SAMPLING  INCREMENT  (Ax)  AND 

HIGHEST  FREQUENCY  OF  INTEREST  (f  ) 

m 
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NUMBER  OF  SAMPLES 


FIGURE  D-2.  NUMBER  OF  SAMPLES  (n  ) VERSUS  LENGTH  OF 

s 

MEASUREMENT  ON  TAPE  (F-F') 
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samples  times  the  sampling  interval.  This  can  be  represented  as: 

^m 

" = (n^-l)'(A-x) 

where 

s = Segment  ratio 

= Length  of  the  time  segment  of  interest  (sec) 
n^  = Number  of  samples 
Ax  = Sampling  increment  (sec) 

Refer  to  Figure  D-3  for  values  of  segment  ratio. 

D,3  ESTIMATE  OF  COMPUTER  TIME  AND  LABOR 

When  the  parameter  values  discussed  in  Section  D.2  have  been  calcu- 
lated and  the  measurements  in  each  group  have  been  further  subdivided  into 
classes  containing  measurements  with  the  same  parameter  values,  an  estimate 
for  performing  a processing  procedure  on  a class  of  measurements  can  be  made, 
and,  from  that,  an  estimate  for  processing  an  entire  group  of  measurements. 

This  is  accomplished  by  selecting  a procedure  to  be  performed.  The 
following  is  a description  of  the  processing  procedures  for  which  guidelines 
have  been  developed  and  a discussion  of  the  method  for  computing  the  corre- 
sponding computer  and  labor  estimate. 

D.3.1  FILTERING  AND  DECIMATION 

This  procedure  consists  of  the  following  functions: 

a.  Retrieve  measurement  file. 

b.  Filter  and  decimate  entire  file  with  respect  to  specified 
frequency  limits. 

c.  Plot  filtered  and  decimated  file. 

d.  Return  entire  filtered  and  decimated  file  to  Archive. 

e.  Retrieve  filtered  and  decimated  file. 

f.  Plot  filtered  and  decimated  file  over  specified  time  limits. 
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SEGMENT  RATIO 


2 3 i*  5 

SAMPLED  PERIOD,  T,  SEC 
T=  (n^-1)  (ax) 


FIGURE  D-3.  SEGMENT  RATIO  (s)  VERSUS  NUMBER  OF  SAMPLES  (n  ), 

SAMPLING  INCREMENT  (Ax) , AND  LENGTH  OF  TIME  SEGMENT  OF 

INTEREST  (T  ) 
m 
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The  algorithm  for  determining  the  uncorrected  amount  of  computer 
time  for  performing  this  procedure  on  a measurement  in  a specific  class  can 
be  expressed  as: 

■ 0.0168  + ^-A_|o.620  +d(0.776  + 0.527  3)1 

where 

t^  = Uncorrected  computer  time  per  measurement  in  a class  (hr) 
n^  = Numbers  of  samples 

d = Decimation  ratio 

s = Segment  ratio 

The  next  step  is  to  calculate  the  uncorrected  amount  of  computer 
time  to  process  a class  of  measurements.  This  is  equivalent  to  the  product 
of  number  of  measurements  in  a class  and  the  uncorrected  computer  time  per 
measurement  in  the  respective  class.  This  can  be  represented  as: 

t = m t 

c cm 

where 

t = Uncorrected  computer  time  per  class  (hr) 
c 

m = Number  of  measurements  in  the  class 
c 

t^  = Uncorrected  computer  time  per  measurement  (hr) 

The  uncorrected  amount  of  computer  time  to  process  a group  of 
files  can  be  expressed  as  the  sum  of  all  the  uncorrected  computer  times  for 
each  class  of  that  group.  This  can  be  expressed  as: 

t =t  +t  +t  

g C-  Cry  C« 


Note:  The  uncorrected  time  will  be  adjusted  by  an  overhead  factor  in  a sub- 
sequent equation. 
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where 


t 

g 

= Uncorrected  computer  time  to  process  a group  of 

measurements  (hr) 

t 

^1 

= Uncorrected  computer  time  for  first  class  of 

measurements  (hr) 

rt 

O 

N) 

= Uncorrected  computer  time  for  second  class  of 

measurements  (hr) 

t 

"3 

= Uncorrected  computer  time  for  third  class  of 

measurements  (hr) 

Finally,  the  computer  time  estimate  for  filtering  and  decimating 
all  the  measurements  is  equivalent  to  0.131  plus  the  sum  of  all  the 
uncorrected  computer  times  to  process  each  group  of  measurements.  This 
can  be  represented  by: 


t = 

0.  131  + t + t + t 

91  9^  93 

where 


t 

Computer  time  to  process  all  measurements  (hr) 

t 

9l 

= Uncorrected  computer  time  to  process  first  group  of 
measurements  (hr) 

t 

^2 

= Uncorrected  computer  time  to  process  second  group  of 
measurements  (hr) 

t 

"3 

= Uncorrected  computer  time  to  process  third  group  of 

measuremen  ts 

DH3 


The  amount  of  estimated  labor-hours  to  perform  this  procedure  is 
equivalent  to  fifteen  (15)  plus  the  product  of  three  and  one-half  (3.5) 
and  the  total  number  of  involved  measurements.  This  can  be  expressed  as: 

h = 15  + 3.5 

m 

where 

h = 1 a bo r- hours  (hr) 

n = Total  number  of  measurements 
m 

An  example  of  the  use  of  this  estimating  guideline  as  applied  to  the 
data  processing  request  described  in  Section  2.1  can  be  found  in  Table  P-3* 

D.3.2  DETREND 

This  procedure  consists  of  the  following  functions: 

a.  Retrieve  measurement  file. 

b.  Filter  and  decimate  entire  file  with  respect  to  specified 
frequency  limits. 

c.  Calculate  integrations  for  motion  measurements  (single 
integration  for  velocities,  double  integration  for  accelerations 

d.  Plot  the  filtered  and  decimated  file  and  respective  integrations 

e.  Return  entire  filtered  and  decimated  file  to  Archive. 

f.  Retrieve  filtered  and  decimated  file. 

g.  Detrend  filtered  and  decimated  file. 

h.  Calculate  integrations  (single  integration  for  air  blasts  and 
velocities,  double  integrations  for  accelerations). 
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I.  Plot  resulting  detrended  file  and  integrations  over  specified 
1 i mi ts . 

j.  Return  detrended  file  to  Archive. 

The  algorithms  for  estimating  the  amount  of  computer  time  for  per- 
forming this  procedure  are: 

(ai r blasts  only) 

'm  - + |0.620  + dd.SH-  l.'.6  s)l 

(velocities  only) 

'm  - ^>1 

(acce lerat ions  only) 

'.n  ■ '>-<>353  + -f;g^S^lo.620  + d(2.87 +1.98  3)1 

t = m t 

c cm 

t =t  +t  +t  

g 

t = 0.l80+t  +t  +t  

9l  % 93 

The  amount  of  estimated  labor-hours  to  perform  this  procedure  is 
equivalent  to  fifteen  (15)  plus  the  product  of  three  and  one-half  (3.5)  and 
the  total  number  of  involved  measurements.  This  can  be  expressed  as: 

h = 15  + 3.5  n 

m 

whe  re 

h = Labor-hours  (hr) 

= Total  number  of  measurements 


D-16 


An  example  of  the  use  of  the  estimating  guideline  for  the  detrend 
procedure  as  applied  to  the  data  processing  request  described  in  Section  2.1 
can  be  found  in  Table  D-A. 


D.3.3  SHOCK  SPECTRA 

This  procedure  consists  of  the  following  functions: 

a.  Retrieve  detrended  measurement  file. 

b.  Calculate  shock  spectra  over  specified  frequency  limits. 

c.  Plot  shock  spectra  on  standard  tripartite  size  grids. 


The  algorithms  for  estimating  the  amount  of  computer  time  for 
performing  this  procedure  are: 


t 


m 


0.00737  + 1.02^4 


n sd 
s 


1 ,000,000 


t = m t 
c cm 


t 

g 


= t 


4-  t 


+ t 


t = 0.0308  + t + t + t 

91  92  93 


The  estimated  number  of  labor-hours  to  perform  this  procedure  is 
equivalent  to  three  (3)  times  the  total  number  of  measurements  to  be 
processed.  This  can  be  expressed  as: 


h 


whe  re 


3.0  n 

m 


h = labor-hours  (hr) 

n = Total  number  of  measurements 
m 


Refer  to  Table  D-5  for  example. 
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TABLE  D-A.  SUMMARY  OF  ESTIMATE  OF  COMPUTER  TIME  AND  LABOR-HOURS  FOR 
DETRENDING  OF  RDA  REQUESTED  DATA 


Gro  up 

t , hr 

g 

Number 

of 

Classes 

Class 

m 

c 

t , hr 
m 

t , hr 
c 

1 

0.1938 

1 

lA 

2 

0.0969 

0.1938 

2 

0.0975 

1 

2A 

1 

0.0975 

0.0975 

3 

0.2541 

2 

3A 

1 

0. 1022 

0.1022 

3B 

1 

0.1519 

0.1519 

A 

0.4530 

1 

4A 

3 

0.1510 

0.4530 

5 

0. 1632 

1 

5A 

1 

0.1632 

0. 1632 

6 

0.1378 

1 

6A 

1 

0.1378 

0.1378 

Total 

Correct 

Time,  t 

= 1.479 

h = 50  Lat 

>or-Hours 

1 
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D.3.^  FOURIER  TRANSFORM 

This  procedure  consists  of  the  following  functions: 

a.  Retrieve  filtered  and  decimated  or  detrended  measurement  file. 

b.  Calculate  Fourier  transform  real  and  imaginary  components. 

c.  Convert  to  amplitude  and  phase  over  specified  frequency  limits. 

d.  Plot  Fourier  amplitude  on  log-log  grid. 


The  algorithms  for  estimating  the  amount  of  computer  time  for 
performing  this  procedure  are: 


tpi  0-00897  + ^ ^000,000 


l4.0  + 0.156  d - 2.06  log^Q(s  d n^) 


t = m^t 
c Cm 


t = t + t + t 
g C3 


t = 0.0377  + t + t + t 

^2  - 


The  estimated  number  of  labor-hours  to  perform  this  procedure  is 
equivalent  to  three  (3)  times  the  total  number  of  measurements  to  be 
processed.  This  can  be  expressed  as: 

h = 3.0  n 

m 


where 

h = Labor- hou rs  (hr) 

n = Total  number  of  measurements 

m 

Refer  to  Table  D-6  for  example. 
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D.3.5  FOURIER  TRANSFORM  AND  SHOCK  SPECTRA 

This  procedure  consists  of  the  following  functions: 

a.  Retrieve  detrended  measurement  file. 

b.  Calculate  Fourier  transform  real  and  imaginary  components 
of  measurement  file. 

c.  Convert  Fourier  transform  to  amplitude  and  phase  over  specified 
frequency  limits. 


d.  Plot  Fourier  amplitude  on  log-log  grid. 

e.  Calculate  shock  spectra  of  measurement  file  over  specified 
frequency  limits. 

f.  Plot  shock  spectra  on  standard  tripartite  size  grid. 


The  algorithms  for  estimating  the  amount  of  computer  time  for 
performing  this  procedure  are: 

n sd  p 

T, "000, 000  + 0.196  d - 2.06  log^Q(s  d n^) 


t = 0.0149  + 


t = m t 
c cm 


t = t + t + t 

g C2  C3 


t = 0.0377  + t + t + t 

9l  92  g 


3 


The  estimated  number  of  labor-hours  to  perform  this  procedure  is 
equivalent  to  three  (3)  times  the  total  number  of  measurements  to  be 
processed.  This  can  be  expressed  as: 


h = 3.0  n 
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where 

h = Labor-hours  (hr) 

n = Total  number  of  measurements 

m 

Refer  to  Table  D“7  for  example. 
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TABLE  D.7.  SUMMARY  OF  ESTIMATE  OF  COMPUTER  TIME  AND  LABOR-HOURS  FOR  FOURIER 
TRANSFORM  AND  SHOCK  SPECTRA  OF  RDA  REQUESTED  DATA 
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